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fluid is calculated for two different theories, Eckart and Israel-Stewart and a comparison 
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velocity solutions that are decreasing exponentially for particular choices of the initial 
conditions. 
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1. Introduction 

Relativistic fluid dynamics provide a simple but intuitive description of many physical 
systems. The hot and compressed nuclear matter behaves like a compressible fluid and 
fluid dynamical effects are observed in high energy heavy ion reactions |]J. Supernova 
or neutron star may be another situation where the relativistic heat-conducting fluid 
gives a satisfactory description. 

The standard theory of the relativistic dissipative fluid proposed by Eckart [[| 
and a similar one by Landau and Lifshitz || was analysed by Hiscock and Lindblom 
@; H; H? 01 an d it was found to be unstable, acausal and ill posed in the linear regime near 
equilibrium. A more complicated extended hydrodynamic theory has been proposed by 
Israel and Stewart || || [L(| which is free of the most part of the troubles of the standard 
theory. 

The Hiscock and Lindblom studies were dedicated to causality, stability, 
hyperbolicity and to the connection between them. In the solution discussion was 
restricted to spatially homogeneous fluid. This paper will improve the dimensional 
restriction considering the evolution equations for the plane symmetric motion (1+1 
dimensional) of the relativistic heat-conducting fluid and propose a complementary 
analysis dedicated to group-invariant solutions study. 

Symmetry analysis is a systematic and accurate way to obtain solutions of 
differential equations. The Lie symmetry group method is well known but the major 
obstacle in the application of this method to partially differential equations is the large 
number of tedious calculations usually involved. 

We begin our analysis by the calculation of the Lie symmetry group for both Eckart 
and Israel-Stewart fluids using, for the equation of state, the high temperature limit of 
the ideal gas, where p = p/3 = nkT. This constraint allows us to eliminate from the 
fluid equations two dependent variables, in particular, we substitute the temperature T 
and the pressure p in favour of the energy density p and the number density n. Next, 
we make a systematic classification of the solutions constructing the optimal system of 
transformations. The investigation of the solutions is completed by solving the reduced 
system of equations; we present both analytical and numerical solutions. 

Section |2| will display the system of equations of the relativistic dissipative fluid 
in its general form eqs. (|l|) and in the particular case of heat-conducting theory eqs. 
(|). The application of the Lie symmetry group method to differential equations and the 
corresponding Lie algebra of the fluid equations (|) are presented in Section ||. Section f| 
is devoted to integrability conditions and invariants, Section |5| to the classification of 
the solutions and Section || to group invariant solutions analysis. Possible applications 
of the group-invariant solutions to high energy collisions, miscellaneous comments and 
final conclusions are discussed in Section |^. 
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2. Dissipative relativistic fluid mechanics 

In the particle frame, used by Eckart, the energy tensor T Q/3 and the particle flux vector 
N a take the form |J: 

T a/3 = pu, a vP +(p + a) (g al3 + u a u p ) + q a u p + q p u a + a a/3 

N a = nu a 

where u a is a unit time-like vector field which may be thought of as the four-velocity of 
the fluid, g af3 is the Minkowski metric diag(-l, 1,1,1), p the energy density of the fluid, 
p the thermodynamical pressure, n the number density, q a is the heat flow, a the bulk 
stress and a a ^ the shear stress. 

The heat conducting case is obtained by setting a = and a al3 = 0. In this 
particular condition the heat flow obeys the following phenomenological law |§: 

q a = -kT (g a ? + u a uP) (t~%T + u^d^up + fau^qp) (2) 

where beside the usual thermal conductivity k and temperature T we have a new 
phenomenological coefficient fa [§]. Setting fa = we obtain the standard Eckart 
heat conducting fluid theory. 
The conservation laws: 

d a T afi = 0, d a N a = (3) 

and the equation (Q) form the complete system of equations of the general theory of 
relativistic dissipative fluids. 

If we write the four- velocity in the form u a = (cosh sinh 0, 0), the heat flow 
vector, which is orthogonal to u a , will be q a = q (sinh cosh 0, 0) and the system of 
equations of the fluid becomes: 

sinh ^>d x n — cosh ^>d t n + n cosh ^d x ^ — n sinh ^d t ^ = 
cosh ^>dtp — sinh ^>d x p + sinh ^>dt q — cosh ^*d x q 

+ [ip + P) sinn ^ + 2g cosh dt$ - [{p + p) cosh * + 2q sinh \&] d x ^ = 

cosh ^>d x p — sinh ^>dtp + sinh ^d x q — cosh ^dtq (4) 

- [{p + p) cosh * + 2g sinh \1>] + [(p + p) sinh \1> + 2q cosh \1>] d x ^ = 

T" 1 (cosh Wd x T - sinh ^d t T) + fa (sinh - cosh ^d t q) 

+ sinh \I/<9 X \I/ - cosh ^/d t ^ + -§f = 

We will supplement the eqs. (f|) by an equation of state (EOS) for the fluid and a 
thermodynamic expression for fa. Choosing the high temperature limit of an ideal gas, 
the equation of state is 

p = p/3 = nkT (5) 

Finally, we take the thermal conductivity to be constant and the second-order 
thermodynamic coefficient fa to be given by fa = 5A/4p; the parameter A will be 1 
for the Israel- Stewart theory and for the Eckart theory. 

The equation (f|) allows us to eliminate T and p in favour of p and n, obtaining 
finally a closed set of evolution equations for four variables (^, n, p, q). 
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3. Symmetry group of transformations and its Lie algebra 

The symmetry group of a system of differential equations is the largest local group of 
transformations acting on the independent and dependent variables of the system with 
the property that it transforms solutions of the system to other solutions. Let 5 be a 
system of differential equations. A symmetry-group of the system S is a local group 
of transformations Q acting on an open subset M. of the space of independent and 
dependent variables of the system with the property that whenever u=f(x) is a solution 
of S, and whenever g ■ f is defined for g G Q, then u = g ■ f(x) is also a solution of the 
system. 

The symmetry group infinitesimal generator is defined by : 

V = rd t + £d x + $<9* + m n + Yd p + VL3 q (6) 

and the first order prolongation of V is: 

pr^V = £d x + rd t + $0* + E<9„ + Td p + Qd q 

+$*<9^ + + Z x d nx + T*d nt (7) 

+v*d Px + v% t + n*d qx + n% t 

where, for example, 

= D x ($ - 01'., - r* t ) + & xx + 7M/„ (8) 

and D x § = $ x + §^ x + § n n x + Q p p x + Q q q x is the total derivative and $ x = d x Q, etc. 

Suppose A v (x,u^ n ') — 0, v — 1,...,/, is a system of differential equations, where 
u ( n ) = (\]> ; n , p, q, ^ x , q tt ). If Q is a local group of transformations acting on Ai and 

p r ( n )y A„(x, w^) = 0, v — 1,...,/, whenever A(x, u^) = 0, for every infinitesimal 
generator V of Q, then ^ is a symmetry group of the system. 

The standard procedure^ is based on finding the infinitesimal coefficient functions 
£, r, $, E,T and fi. Substituting the general formulae for <£> x , T, x , etc. and equating 
the coefficients of various monomials in the first and second order partial derivatives 
of ty,n,p and q, we find the defining equations. We wish to determine all possible 
coefficient functions £, r, E, T and fl by solving the defining equations system so that 
the corresponding one-parameter group exp(eV) will be the symmetry group of the 
equations (f|). These solutions are: 

Eckart : r = c\ + tc^, £ = c 2 + XC4, $ = 0, E = — nc^, T = pc%, Q = qc% 

(9) 

Israel — Stewart : r = c\ + £04, £ = c 2 + XC4, $ = 0, E = 0, T = pc 3 , f2 = qc% 
where constants. The basis of the corresponding Lie algebra are: 

Eckart : V x = d t , V 2 = d x , V 3 = pd p + qd q , V 4 = td t + xd x - nd n 
Israel — Stewart : V\ = d t , V 2 = d x , V3 = pd p + qd q , 

I A good description can be found in Q 
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4. Solvable group and invariants 

Because we have the Lie algebra of the equations (^) we want to know if the general 
solution of the system of differential equations can be found by quadratures. This thing 
is possible if the Lie group is solvable. The requirement for solvability is equivalent to 
the existence of a basis {Vi, . . . , V r } of Lie algebra g such that 

[^] = E4 y * ( n ) 

k=l 

whenever i < j. Looking at the commutator table of the Lie algebra we will see that 
the requirement of solvability is fulfilled for both Eckart and Israel-Stewart theories. 

Table 1. Commutator table for the Eckart and Israel-Stewart algebra 
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We use the method of characteristics to compute the invariants of the Lie algebra 
hopping that the reduced system, which can be obtained using the invariants of the 
group, will help us to solve the system of equations (|J). An n-th order differential 
invariant of a group G is a smooth function depending on the independent and 
dependent variables and their derivatives, invariant on the action of the corresponding 
n-th prolongation of G Suppose that we have the following generator: 



V = nd t + £id x + $ i( 9* + Ed n + T.dp + Qid g (12) 

A local invariant ( of Vi is a solution of the linear, homogeneous first order partial 
differential equation: 

V(0 = Tid t ( + + $i9»C + EAC + FidpC + CW = (13) 

The classical theory of such equations shows that the general solution of equation 
(|13|) can be found by integrating the corresponding characteristic system of differential 
equations, which is 

dt dx dty dn dp dq 

n o 'i', Si r, Qi 

Doing this integration we get, in this case, five invariants; we now re-express the next 
generator of Lie algebra in terms of these five invariants and then we perform another 
integration. We continue this calculation until we re-express and integrate the last 
generator; at this point we obtain a set of invariants that represent the system of 
independent invariants of this group. The system of invariants can be used to reduce 
the order of the original equations - constructing the reduced order system of equations. 
Doing this one can hope to find simple equations that can be integrated (for example 

H). 
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5. Classification of group-invariant solutions 



A solution of the system of partial differential equations is said to be (/-invariant if it 
is unchanged by all the group transformations in Q. In general, to each s-parameter 
subgroup Ti of the full symmetry group Q of a system of differential equations, there 
will correspond a family of group-invariant solutions. Since there are almost always an 
infinite number of such subgroups, it is not usually feasible to list all possible group- 
invariant solutions of the system. We need an effective systematic means of classifying 
these solutions, leading to an optimal system of group-invariant solutions from which 
every other solution can be derived. Since the elements g G Q not in the subgroup Ti 
will transform an 7i-invariant solution to some other group-invariant solution, only those 
solutions not so related need to be listed in our optimal system. An optimal system of 
s-parameter subgroups is a list of conjugancy inequivalent s-parameter subgroups with 
the property that any other subgroup is conjugate to precisely one subgroup in the list 
(conjugancy map: h — > ghg~l). Let Q be a Lie group with Lie algebra g; for each v G g, 
the adjoint vector ad v at w G g is ad v\ w = [w, v] = — [v, w]. If v G g generates the one- 
parameter subgroup Ti = [exp(ev) : e G 9ft], then Ad g{v) will generate the conjugate 
one-parameter subgroup gTig^ 1 . Now we can construct the adjoint representation Ad Q 
of the Lie group by summing the Lie series 

Ad(exj>(ev))w 



—Aad v) n (w) =w- e[v, w] + ^—[v, [v, w)) 

n=0 n - 1 

obtaining the adjoint table. 

Table 2. Adjoint table for Eckart and Israel-Stewart algebra. 
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Finding the optimal system is a well known standard method JTT|] ; we have found 
an optimal system of one-dimensional (a G 3ft) sub-algebras to be that spanned by: 



Eckart 

d x 
d t 

xd x + td t - nd n 
d t + ad x 

xd x + td t — n d n + a pd p + aqd q 
pd p + qdg + ad x + d t 
pd p + qdg + d x 
pdp + qd q + d t 
pdp + qdg 



Israel — Stewart 

d x 

d t 

xd x + td t - nd n 
d t + ad x 



(16) 
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6. Group-invariant solutions study 

Each symmetry of the optimal system of transformations (|T6| ) will be treated separately. 

1) The fluid is spatially homogeneous and because the reduced system is very 
simple the fluid equations can be immediately integrated analytically The Figure [I] 
shows that, for any initial state, the Eckart theory is unstable: as t — > oo then 
v — > \ => T — > oo, n — > 0, p = 3p ^ 3E and q — > —2E . The same thing happens, 
for the Israel-Stewart theory, only if the initial velocity is higher than the critical value 
0.76 but it can be stable for initial values v < 0.76, where v = tanh\l/. This spatially 
homogeneous case was also studied by Hiscock and Lindblom [|7| and they found the 
critical value around v c = 0.51. 




Figure 1. Israel-Stewart and Eckart velocity solution, where t is ' ° t. 

2) Also in this stationary case, the reduced system is very easy to integrate 
and the space dependence of the velocity can be observed in the Figure |2|. For 
the Eckart fluid, if the initial velocity is higher then 0.77 then as x —>■ oo we have 
v — >l=^T^oo, n — > 0, p = 3p — > 3E , q — > —2E . If we correlate this stationary 
behaviour with the previous result, we can comment that the stationary solution will 
increase in time for any initial Eckart velocity conditions; the Stewart-Israel fluid has 
the chance to decrease when t — > oo, x — > oo for a initial velocity smaller that 0.76. 




0.60.8 1 1.21.4 2 4 6 SIO 



Figure 2. Israel-Stewart and Eckart velocity solution, where x is ^x. 

3) The invariants of the characteristic equation dx/x = dt/t = —dn/n are 
a = n ■ t, y = x/t and for the Eckart fluid the reduced system of equations take 
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a 1 a 



„ — — j- [In (sinh ^ + y cosh \E0] 
ay 



a 



N (sinh * + y cosh ^) 1 



p y (1 + y 2 ) cosh(2^) + 2(-l + y 2 + y sinh(2^) 



6*, 



1-f 



2ycosh(2fr) + (l + ?/ 2 )smh(2fr) 
+ 4p* ?/ + p y : = 



;i7) 



3fciV n 



Py 



i-y 2 

1 + [(1 + y 2 ) cosh(2^) + 2y sinh(2^)] 







X sinh \1> + y cosh \1> " r (1 + y 2 ) sinh(2\l>) + 2ycosh(2vI/) 

where = ^,etc. The system of equations can be decoupled and the behaviour of 
the velocity numerical solution can be observed in the first picture of the Figure EjL 



Eckart 




Eckart 




Figure 3. Eckart velocity solution for the 3 and 4 optimal transformations 



It is important to mention that the first picture of the Figure |3| shows that we have 
both exponentially increasing and decreasing solutions, which means that we found a 
good Eckart velocity solution. Also important is that if the initial velocity is smaller 
than 0.75 we have only decreasing solutions. Further numerical investigations of the 
equation are needed in order to obtain information about the stability, initial condition 
dependence, etc.. 

For the Israel-Stewart fluid the reduced system is much more complicated and we 
didn't succeed to solve the equations. 

4) For this optimal transformation the invariant is y = t — ax; if we choose a = — 1 
the reduced system can be easily decoupled and for the Eckart fluid the solutions are: 



tanh 



In 



v / X Ci/(fciV )tanh U X Ci/(kN )(y - C 2 



n = N oy J(l + v)/(l-v), p y = 0, q= -^-p 



;is) 



3kN r (l + v) 2 

The velocity solution is ploted in the second picture of the Figure [| This t + x 
dependence is typical for a wave type solution but asy = t + x— >oowe arrive again to 
wrong physical situation: v —>■ 1 =^ n — > oo, T^O. The choice a = — 1 does not affect 
the physics, it can change the x axes direction or it can scale the space coordinate. 

The Israel-Stewart reduced system is also integrable but we found complex velocity 
solution. 
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5) The invariants are (3 = n-x, y = x/t, w = p-t a , 6 = q/p. After some calculations, 
the reduced system of equations becomes: 

(3 = N y/(siah * + y cosh *) 

= -{(In [(sinh * + y cos h *) (cosh * + y sinh )]) y + + ysinh vl> 

+ a((l + y) sinh(2*) + 2y) 

(sinh* + y cosh*) (cosh* + y sinh *)(6(y 2 - 1)*^ - 3a) 

| 3JfeJV 2ycosh(2*) + (1 + y 2 )cosh(2*) + 2(y 2 - 1) j _1 ( 19 ) 

X (sinh* + y cosh*) (cosh* + y sinh*) (6 (y 2 - l)* y - 3a) 
_ [w^Wy (2y cosh 2* + (1 + y 2 ) cosh 2* + 2(y 2 - 1)) + a ((1 + y) sinh 2* + 2j/)] 
~ (6(y 2 -l)*,-3a) 

-2 (y sinh 2* + cosh 2* - 2 - 3yd) - 3y6 y + a (sinh 2* -36) -2 (2y + 36) * y = 

The system of equations can be decoupled and solved numerically. The behaviour of 
the numerical velocity solution is the same as that was founded for the third optimal 
transformation; ofcourse, more precise numerical investigations are also needed. 

6) In this case the invariants are y = x ± at, a = pexp (—t), a = q/p and the 
reduced system is: 

n y (sinh * — a cosh *) + n (cosh * — a sinh *) * y =0 
o~ x o- y [a cosh * — sinh * + aa sinh * — a cosh *) + a y (a sinh * — cosh *) 
+ cosh * + a sinh * + [4/3 (a sinh * - cosh *) + 2a (a cosh * - sinh *)] * J/ = 
o~ x o y (1/3 cosh * — l/3a sinh * + a sinh * — aa cosh *) 
+a y (sinh * — a cosh *) — 1/3 sinh * — a cosh * 
+ [4/3 (sinh * - a cosh *) + 2a (cosh * - a sinh *)] * y = 
cosh * (p~ 1 pa~ 1 a — n~ 1 n y ) — sinh * [p~ 1 p (acr^V + 1) — n~ l n y \ 

+ (sinh * - a cosh *) * y + 3nak/n = 

We can decouple the equations and in the end we obtain the following equation: 

* ra + k^ 2 y + k 2 f (*) *, + k 3 g (*) = (21) 

where ^1,2,3 are constants and / (*) , g (*) are known functions of *. The equation (pl|) 
is integrable [12]; because the integration algorithm is very long we will present only the 



final form of the equation that can be integrated: 

d*(y)/dy=-t£'{t)exp(C 2 V) (22) 

where £ = C^exp [(2 — C2)*] and Cj are constants; t, which is a function of * via £, 
has to be found from the equation t 2 ^"(t) +C^(t) = 0. Because the analytical solution 
formula is very long and complicated we don't write it here, we prefer to mention that 
it has the same exponentially increasing trouble founded before. 

For the 7), 8) and 9) transformation we can't reduce the fluid system of equations. 
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7. Summary and outlook 

The Lie symmetry group method was applied to two different relativistic heat- 
conducting fluid theories, the Eckart theory and the Israel-Stewart theory, and a 
comparison between them has been made. We may summarise the results as follows: 

• applying the Lie symmetry group method to both fluids, we found that the Lie 
algebra of the Eckart theory has an additional transformation, namely the space, 
time and particle density dilatation (V4). 

• a systematic classification of the solutions is done by constructing the optimal 
system of transformations. 

• solving all the Eckart 's reduced systems of equations, we found that during the 
evolution of the fluid, there are, beside the exponentially increasing velocities, also 
solutions with good physical behaviour; we mention that in the 1 dimensional 
approximation, for the Eckart theory, there are only exponentially increasing 
solutions. 

• due to a different (more complicate) phenomenological law, the optimal 
transformation system of the Israel-Stewart theory is trivial and as a consequence 
the reduced systems of equations are not very much simplified but the 
transformation invariants; anyway, we found good physical behaviour for the 
spatially homogeneous and stationary solutions. 

Due to the fact that each fluid has a different phenomenological law, the 
corresponding symmetry groups is different. We considered here only two values for 
the parameter A that leads us to some well known fluids, but performing the same 
kind of analysis for other A values, we could shed some light on the domain in which 
the relativistic heat-conducting fluid theory is well behaved. Taking into account the 
bulk stress and the shear stress, the phenomenological law becomes more general but 
then the fluid equations will be more difficult to solve and in spite of the tremendous 
work needed, it is possible to obtain a very simple symmetry group that does not help 
very much on the reduced system integration. It is important to point out that we are 
looking to find the simplest well behaved fluid theory in order to be solved much easier. 
Anyway, it is an open subject that needs some attention. 

Because the 1+1 dimensional approximation is appropriate to describe the collision 
of two highly Lorentz contracted heavy ions, further analysis of the relativistic fluid, 
with or without viscosities, with different A values, etc., are very important because, for 
some years, any theoretical predictions can be compared with the experimental data that 
were taken at AGS or CERN. The concept of a quark-gluon plasma (QGP) predicted for 
heavy ion collisions [H| has been the driving force for the experimental and theoretical 



studies of hadronic matter at high energy densities. A lot of theoretical work has been 
done in order to use the hydrodynamical models in heavy ion physics and several models 
based on relativistic hydrodynamics exist, the reader is reffered to Refs. 0, 14, 15]. 
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